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Using a generalized T-matrix description which, in principle, exactly includes Coulomb correla- 
tions and potential scattering events, resonant and bound impurity states are discussed. Like in the 
non- interacting case, the effects of the scattering potential can be divided into different partial wave 
channels, exploiting the symmetry of the underlying lattice. Due to Coulomb correlations bare local 
(i.e. s-wave) potentials become dynamic and extended, being responsible also for p-, d-wave etc. 
scattering effects. Numerically exact results for both the two-dimensional t-J and Hubbard mod- 
els are used to construct a simple (static) approximation to the effective impurity potential which 
is shown to reproduce the exact resonant scattering and bound states in the relevant symmetry 
channels. 



Experimental data on high-Tc materials including sub- 
stitutional defects serve for probing their unusual nor- 
mal state properties as well as the nature of the super- 
conducting phase. Introducing, for example, Zn in the 
copper oxide planes influences transport and magnetic 
properties and leads-|tp| a drastic reduction of the transi- 
tion temperature Tjilcl. The physical conclusions drawn 
from the results are still controversial. For example, this 
is the case for the local magnetic moment induced by a 
non-magnetic imputit}j-|Such as Zn in hole doped super- 
conducting samplealla wbich could lead to a magnetic 
pair breaking mechanisrrC'Q. 

However, one common feature seems to be the spa- 
tially extended nature of the magnetic effects on the 
CuO-based high-Tc compounds. From residual resis- 
tivity measurements one deduces large scattering cross 
sections of several lattice constants in diameteiaa favor- 
ing strong potential scattering likely connected to d-wave 
pair breakingH. On the theoretical side there has been 
work based on s-wave impurity scattering in anisotropic 
superconductors revealiag the strong influence of non- 
magnetic impuritiest2rtil. In addition, calculations were 
carried out for d-wave superconductors with conabippd 
scattering by impurities and Coulomb correlationala'ta. 

In this paper we focus on the interplay between strong 
correlations and local potential scattering. This situation 
can be ascribed to Zn^+ impurities with a filled (Sd)^*^ 
shell in the cupratc superconductors. We derive a micro- 
scopic picture of the mediated scattering processes using 
a generalized T-matrix description which exactly includes 
Coulomb correlations and which exploits the symmetries 
of the underlying lattice. Partial-wave phase shifts are 
introduced which, like in the non-interacting case deter- 
mine the additional density of states caused by the im- 
purity potential. For this purpose we take advantage of 
numerical exact results of translational-invariant lattices 
and combine them in analytical T-matrix expressions for 
a single impurity. By doing this we can determine the 



effective potential strengths, acting in a perturbed t-J 
system investigated earlierEZl 

Consider a correlated two-dimensional system such as 
the one band Hubbard model with Hamiltonian, 



H = Hn 



U 



^n,;|n,x, (1) 



using standard notations, and a bare impurity operator 



Vi. 



7np,a 



> ybare 



(2) 



Here the four bonds connected to the impurity site at 
the origin are modified by toa- The simplest description 
of a vacancy is then obtained by setting ioo- = t- For 
the sake of simplicity the potential is restricted here to 
the central site but the following analysis is more general 
and valid for extended one-body potentials. The above 
scattering problem on a lattice is defined spin-dependent 
and can therefore be adapted to different physical scenar- 
ios. In the case of non-interacting systems, e.g. U = 0, 
successive scattering is described by the Greens function 
Ga = Go + Go Ta Gq. Here, the T-matrix, defined as in 
usual scattering theory, T„ = Vimp,a (1 - Go Vimp,a)~^, 
involves only the unperturbed propagator Go and the 
bare potential Vimp,a- 

On the other hand, for correlated systems it is ad hoc not 
clear how to define a T-matrix leading to a similar result 
for the exact Greens function G. This can be achieved 
by investigating the representation of G extracted from 
perturbation theory. For finite temperatures, the exact 
Greens Function G can be written in the Matsubara for- 
malism and in momentum space as 



G-i(k,k',c^„ 



e(k) -E,(k,k',w„ 



(3) 



The self-energy E is a sum of three terms I]c/(k, k, a;„), 
S(7y,cr(k, k', a;„) and V/mp,(T(k, k'). In a diagrammatic 
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analysis S[/(k, k, ujn) contains exclusively many-body ef- 
fects, the mixed self-energy term I](7v_o-(k, k', aj„) in- 
cludes all correlation diagrams-with an internal line scat- 
tered at least once by Vfmp.crtS, examples of which are 
given in Figs, lb) and Ic). V/mp^^lk, k') is then just a 
single scattering event caused by the bare impurity po- 
tential. Defining a dynamical effective potential 

Ve//^<,(k,k',cj„) = Ec/K<T(k,k',w„) + V^/,„p,<,(k,k') (4) 

one arrives at the desired formulation, i.e. 

Gcrlkjk^Wn) = Gc/lk, k,w„) 

+ Gc/ (k, k, w„ ) (k, k' , ) Gt/ (k' , k' , a'„) , (5) 

where G^"'^(k, k, aj„) = iw„ — e(k) — ^[/(k, k, a;„) and, by 
omitting the momentimr indices, 

Tcr(w„) = Veff,cr{(^n) (1 - Gc/(w„) Vef f^a{l^n))~^ ■ (6) 

In order to investigate the one-particle behavior of the 
system by means of the T-matrix, one exploits the trans- 
formation of its matrix elements in real space according 
to the symmetrized states of the underlying point-group 
symmetry of the latticdl^. This means that the effects of 
the scattering potential, which respects this symmetry, 
can be divided into contributions to different partial- 
wave channels. Note that for the non-interacting case, 
the bare potentials Vq^^"^ and toa only contribute to s- 
wave processes, while scattering in the d-, p-wave etc. 
channels would require extended bare potentials. As seen 
later, this is substantially modified by the Coulomb in- 
teraction leading to correlation induced longer-ranged ef- 
fective potentials. The poles of the retarded T-matrix, 
which correspond to the zero-values of its determinant, 
either determine resonant states, which are found inside 
the bands of the unperturbed {Vimp = 0) system, or 
bound states separated from the continuum. Xhese state- 
ments are clear for non-interacting systemalj but hold 
also for U ^ 0. Note that Vefj^ai^^) is a real quantity 
at the bound state energies lu ~ lobs due to the infinite 
lifetime of these states. 

By using the irreducible representations the complex 
determinant of the T-matrix factorizes into the subdeter- 
minants of its partial-wave decompositions. The partial- 
wave phase shifts 4>a{^) are then introduced as the (neg- 
ative) phase of these a-wave subdeterminants. They de- 
termine the additional density of [States (DOS) caused 
by an extended impurity potentiaLJ. The total DOS 
of the system is given by the trace over the retarded 
Greens function G'"'^* and can be written as the sum of 
Du{lo), the DOS of the correlated system without the 
impurity and an additional impurity-caused contribution 
AZ?cr(i-^). Bringing in the phase shifts, this leads to 

l^D,{uj) = -^Im^ In Det(l - G^ (^) ^//..(c.)), 
= EAi.„,.(.) = ^5:.,(a)^, (7) 



where A'' counts the number of sites and g{a) denotes the 
degeneracy of the representation a. 

This shows that the phase shifts and the T-matrix 
defined for non-interacting systems are still helpful con- 
cepts for studying bound states and resonant scattering 
in the interacting case. However, the use of Eq.(^) im- 
plies the knowledge of the effective potential which is a 
priori frequency and interaction (U) dependent. Note 

that Veff can be expressed as T4//,ct('^) = G'^''^{lo) — 
G'''J'^{lo) ^, which, in principle, would make it possible 
to extract the dynamical potentials by numerical calcu- 
lations. Here we adopt a somewhat simpler approach by 
showing that a crude static but nevertheless interaction- 
dependent Ansatz can be made for Veff^a which can de- 
scribe, within reasonable accuracy, the scattering pro- 
cesses. This method enables us to understand in more 
physical terms the interplay between the interaction and 
short-range impurity scattering. 

One important feature which has to be included in the 
Ansatz is the longer range character of Vef / due to the 
background of the correlated host. This can simply be 
seen from the connection between the bare and the dy- 
namically modified potential. By explicitly transforming 
the Coulomb term to the real-space irreducible repre- 
sentation, one finds local interactions between different 
wave-symmetries, i.e. 

V^QLx ,a2 ,CK3 ,CK4 

where c"lJ^ are the symmetrized annihilation (creation) 
operators of the one-particle spin (a) states, a labels 
in-general all irreducible representations of the site-type 
J/E3. Note that, although it contains a-wave operators 
the above expression for Hij is invariant under all the 
point group symmetries. This implies selection rules on 
the matrix elements {7^i"2a3a4 -^hich^ for simplicity, are 
not specified here. Therefore, even for a pure s-wave po- 
tential V^"-''' J2a CctqCcto a.t the origin, incoming d-wave 
propagators are scattered at v = 1 from a potential me- 
diated by the Coulomb repulsion as shown in Fig. 1. 

This demonstrates that, due to many-body effects, 
even local (s-wave) potentials always become dynami- 
cally extended, being in general responsible for scattering 
contributions in all symmetry channels. In our numer- 
ical results the effective potential in Eq.(4) is described 
by assuming the following static approximation for the 
effective potential in Eq.(4); 

Generally, the j/-off-diagonal (diagonal) part of Veff 
also takes into account the effective change of the hop- 
ping amplitudes (potential terms). 

Next, we proceed with our numerical results. At zero 
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temperature, the retarded Greens function Gjj appear- 
ing in the above equations is calculated by exact diag- 
onalizalion (ED) techniques for the half-fillecL^D t-J 
modeled. In this context, we regard this modelO as the 
strong-coupling limit of our original Hamiltonian Eq. (0) 
to perform the T~matrix calculationS|-(i.e. J = At^ /U). 
We consider the situation of referenccED, where an inert 
site in a t-J lattice produced bound states of different 
wave symmetries a in the singlet one-particle excitation 
spectrum. The above formalism is used to determine the 
effective potential strengths causing these bound states. 
The description in terms of spin-dependent (effective) 
potentials takes into account that an excess-spin lifts the 
spin-degeneracy of the system, leading also to an antifefc. 
romagnetic structure in the vicinity of the impurityli3'E3 
In contrast, an infinitesimal magnetic coupling of the 
inert site to the lattice restores the spin-degeneracyEj, 
which would lead again to spin-independent effective po- 
tentials. Keeping in mind, that we restrict our calcula- 
tions to the low^ying excitations in the above mentioned 
singlet channeO, we omit spin indices from now on. For 
the case of interest here (t^ = t or more generally for 
ybare _ -(.j^g problem bccomes particle-hole symmet- 
ric. Then, for simplicity, we define a; = as the lower 
edge of the quasiparticle band of the J-model (lower 
Hubbard band). In the hole representation, positive (neg- 
ative) values of V^i 2 correspond to repulsive (attractive) 
potentials. ■— ■ 

According to the local DOS results of Ref O, the values 
of VJf and are set to large values ~ 20i in order to ex- 
pel particles from the corresponding site-types u = Q,2. 
In our calculations t' increases the bound-state energy 
qixajdratically. Taking the actual small binding energies 
oillZl into account the small parameter t' is therefore set 
to zero. In this sense the derived potential strengths rep- 
resent an upper limit. 

Fig. 2 displays the exact bound state energies u}BS^ cal- 
culated from the local DOS for a 20-site cluster with an 
isolated site at the origiiO. By using the T-matrix, the 
energetic locations of the bound states are reproduced 
for values of J = 0.5 and J = 1.0 fixing the effective 
potentials . Fig. 2 also exhibits the dependence of 
on the exchange coupling J of the unperturbed system. 
As one expects, the attractive potentials increase with 
J. One finds the actual value of the effective potential to 
be of the order of J, which is the loss of magnetic energy 
per bond of the inert site. 

Due to the dynamic renormalization, the potential 
acts, as mentioned, differently in each symmetry sec- 
tor being largest for the s-wave channel. Note that in 
the limit of vanishing interaction strength ([/ — > 0), the 
p- and c?-wave potentials also have to vanish, reflect- 
ing non-extended potential scattering. Therefore, in an 
interacting system, there should exist an optimum corre- 
lation strength producing the largest effective potentials 
for bound states. 

For J = 0.5 Fig. 3a) shows the DOS of the unper- 



turbed 26-site cluster used to calculate the additional 
impurity-induced density in Fig. 3b). Assuming the 
effective potentials to be slowly varying functions of uj 
one finds a suppression of denfity for small energies in 
agreement with the exact resultO in Fig. 3c). 

Now, we investigate scattering in the 2D Hubbard 
model. For finite temperatures, we exploit results 
of Quanturp_.Monte Carlo (QMC)/Maximum Entropy 
calculationsEj for the undoped system by inserting the 
unperturbed Greens functions into the T-matrix equa- 
tions. 

At half-filling and for U = 8i, (3t — 10, we use 
the effective potentials derived from the i-J-model for 
J — At'^/U = 0.5 and for the particle-hole symmetric 
case to — t. In this case, Vq'^^^ is irrelevant and can be 
set to zero, while the effective potentials change sign at 
the chemical potential /i. Both pure systems are char- 
acterized by long-range antiferromagnetic correlations 
being disrupted by the impurity which should lead to a 
similar behavior regarding the formation of bound states. 
Actually, the states appearing below the t-J quasiparti- 
cle band are now located just above the small dispersive 
quasiparticle band of width ~ 2 J riding on top of a sev- 
eral t wide incoherent backgrounder. This can be seen 
in Fig. 4 displaying the unperturbed DOS Du obtained 
from QMC calculations and the (non-normalized) ad- 
ditional density AD a- Differing from the i-J-case the 
s-wave contributions are shifted towards the quasiparti- 
cle band. The appearance of these states results in an 
overall effect of narrowing the gap of the unperturbed 
Hubbard model by almost 20%. 

In summary, we have considered a formally exact T- 
matrix method adapted to correlated systems which 
revealed the many-body interactions as the origin of dy- 
namically extended potentials. Using numerical data for 
the unperturbed propagators, we reproduced the static 
effective potentials causing bound states in a 2D t-J- 
lattice with an inert site. Applying the formalism to the 
insulating 2D Hubbard model, we have shown the corre- 
spondence of the two models regarding the appearance of 
spectral weight in the correlation gap. Our result is rele- 
vant for T-matrix approximations in the dilute impurity 
limit, i.e. where the single impurity scattering has to be 
treated exactly. There, higher symmetry channels than 
the usually considered s-wave contributions have to be 
taken into account for interacting systems. Our static, 
effective potentials in Eq.(9) can, for example, directly be 
used as an input to a cluster diagonalization approach, 
where the effect of potential impurity scattering in a 2 11 
superconductor doseribed by a BCS mean-field theorjO 
is studied. In ref.Ej the impurity potential was param- 
eterized. Nevertheless, already this work demonstrated 
that the range of the impurity potential is of quantita- 
tive importance in the case of strong potential scatterers. 
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FIG. 1: Typical lowest order diagrams for 
the self-energy: a) Local interaction ver- 
tex connecting different symmetry channels. 
b),c) Diagrams of T^uv describing (/-wave ex- 
tended potential scattering of type I; first or- 
der static contribution (b) and second order 
dynamic scattering process (c) in U. Scat- 
tering by the bare potential is depicted by a 
cross. 
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EIG. 2: Bound state energies lobs of an inert site in the t-J-modei as a 
function of J [17]. Tlie effective potentiais V-^^ causing tliese states, derived 
by tlie generaiized T-matrix metliod from a 26-site cluster. 




FIG3: a)Total density of states of the 26-site cluster without the impurity 
calculated exactly by ED. b) Additional density due to effective potentials at 
site types i/ = 0, f , 2. d-wave: short-dashed lines, s-wave: solid lines, p-wave: 
long-dashed lines, c) Exact local density of states at site type i/ = f for a 
20-site cluster with an inert site [f7]. 
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FIG. 4: Influence of effective potentials on the half-flUed Hubbard model. 
Du from QMC calculations [(3t = 10). The appearance of various symmetry 
bound states causes gap-narrowing at the edges of the Hubbard bands. 



